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o 

r^ , We present colours of sources detected with the Herschel/SPIRE instrument in deep extragalactic surveys of the Lockman Hole, Spitzer-FLS, 
O !■ and GOODS-N fields in three photometric bands at 250, 350 and 500 pm. We compare these with expectations from the literature and discuss 
associated uncertainties and biases in the SPIRE data. We identify a 500 /jm flux limited selection of sources from the HerMES point source 
catalogue that appears free from neighbouring/blended sources in all three SPIRE bands. We compare the colours with redshift tracks of various 
contemporary models. Based on these spectral templates we show that regions corresponding to specific population types and redshifts can be 
^ . identified better in colour-flux space. The redshift tracks as well as the colour-flux plots imply a majority of detected objects with redshifts at 

1 ^i ' l<z<3.5, somewhat depending on the group of model SEDs used. We also find that a population of S250/S350 < 0.8 at fluxes above 50 mjy as 

observed by SPIRE is not well represented by contemporary models and could consist of a mix of cold and lensed galaxies. 

' Key words. Submillimeter: galaxies - Galaxies: evolution - Galaxies: high-redshift 

fT^ 1 . Introduction tral energy distribution (SED) from these galaxies is redshifted 

Cn| ^ , ,. . , , . • • . , into the SPIRE bands. Deep SPIRE surveys will thus provide a 

. Galaxy formation and evolution are major topics in cosmology ^^^ ^^^^^^ ^^ ^^^ ^^^ -^^^^^ j^^j^^^ especially for 

in and asti-ophysics. Recently, sophistic ated computer simulations ^j^^ ^^^^ energetic dust-rich objects. 

^ based on hierarchical formation (eg. ( |Sprmgel et al. || 2005 | ) ) and j„ ^j^^^ ^^ investigate SPIRE colours in three bands, 

^ deep cosmological surveys across the electroniagnetic spectrum 250 ;,,„, 350 ^/m, and 500 tim, of a sample of sources selected 

. have transformed the field from a collection of a few hypotheti- ^^ 5QQ -^ ^^^^ ^^-^^^^ Demonstration Phase (SDP) fields of 

_ > . cal models to a mature science that is profoundly shaping our ^j^^ ^^^^^j^^j Multi-tired Extragalactic Survey (HerME^ key 

[V : understanding of the univer se. The Spectral a nd Photometric p^ject (Oliver et al. 2010). Our goal is to constrain the evolution 

^ Imaging Receiver (SPIRE, griffin et alJj2gl^, on board of ^^ ^j^^ ^^ ^^^ ^^ FIR/submm galaxies (i.e. the SED vs. redshift 

C^ ESAs Herschel observatory (Pilbratt et alj|2010|), opens a new relation). In particular we want to examine whether there is a 

window between 250 - 500 A/m for observations of heavily dust ^^^ population of submm galaxies whose SEDs do not match 

obscured high redshift (z ~ 2 - 3) galaxies. The combination of ^j^^ predictions of cun-ent models. 

sensitivity and high resolution provided by its 3.5m telescope, 

the largest ever launched, allows large area confusion limited 

far-infrared and submillimeter (FIR/submm) surveys to reach 

unprecedented depths, exploring previously undetectable remote 



o 



2. Observations and data processing 



galaxy populations. Studies of the integrated cosmic background The HerMES key project is constructed in order to obtain a com- 
light have shown that at least half the radiation from all galaxies plete bolometric census of star-formation in the Universe. It con- 
lies in the FIR/submm. In particular, galaxies at z > 1 radiate sists of 6 tiers of survey fields with increasing depth over smaller 
mosdy in this band Lagache et al. (2005). The peak of the spec- areas, covering most of the fields on the sky observed across the 



electromagnetic spectrum by state-of-the-art facilities plus indi- 

* Herschel is an ESA space observatory with science instruments yidual selected clusters. A total of 4 HerMES fields were sur- 

provided by European-led Principal Investigator consortia and with im- 

portant participation from NASA. ' hermes.sussex.ac.uk 
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veyed during Herschel's SDP and we have used the deep ob- 
servations in GOODS-N, Lockman-North, FLS, and Lockman- 
SWIRE for our analysis. The covered areas are 0.25, 0.34, 5.81 
and 13.2 deg^ respectively with relative depths of 1.0, 0.23, 0.05, 
0.033 that were calculated as the fraction of the number of re- 
peats and scan speed, normalised to the deepest fie ld GOODS- 
N. M ore details of the observations are given by lOliver et aP 
dMO). 

Data processing based on the standard SPIRE Scan Map 
Pipeline JGriffin et all l2008h yielded maps in the three SPIRE 
bands, and source catalogues for each indiv idual band were 
gene rated using the SUSS EXtracto r software ( Savage & Oliveij 
12007.) within HIPE 3.0 ( OttetalJ 12006). The three shallower 
maps were smoothed with point-source optimised filters while 
the deepest map was filtered with a delta function to find sources 
and separately with a 3 x 3 pi xel point so urce response func- 
tion (PRE) to extract the fluxes (lOliveret alj201 0). The FLS and 
Lockman-SWIRE fields were Wiener filtered to reduce effects 
by diffuse Cirrus. For the source extraction a Gaussian PRE was 
assumed, with FWHM of 18.2", 25.2" and 36.3" for the SPIRE 
250, 350 & 500 jj m. bands, respectiv ely. Details of the proce- 
dure are given by lOliver et alJ (1201 Ol) and ISmith et alJ (l2010l) . 
They attained formal 1-cr point source uncertainties of 5.7, 7.4 
and 7.8 mJy for GOODS-N, 7.0, 8.5 and 8.8 mJy for Lockman- 
North, 9.0, 10.3 and 10.6 mJy for FLS, and 11.1, 16.9 and 
15.1 mJy for Lockman-SWIRE, respectively. These numbers in- 
clude a contribution from source confusion of approximately 
5.6, 7.4, and 7.7 mJy for GOODS-N, 6.8, 8.3 and 8.5 mJy for 
Lockman-North, 8.4, 9.8, 9.7 mJy for FLS, and 8.5 , 13.9, 10.2 



for Lockman-SWIRE. For the two deepest fields Smi th et al 
(2010) attribute the differences to the results of Nguven et al 
(1201 Ol) . and the differences between the fields mainly to the 



source extraction method used. 

The catalogues contain additional parameters to allow for 
quality checking and source selection. These are: /) A formal 
eiTor in the flux measurement, propagated through from the er- 
ror maps created by the map-maker, representing a fair estimate 
of the instrumental noise. //) A total error that is the quadratic 
co-addition of instrument noise and average estimated confusion 
noise over the map. ///) Two separate flux estimates for the same 
positions using two different halves of the data (half-maps), sep- 
arated in time, allowing for the detection and exclusion of spuri- 
ous sources, mostly due to high energy particle hits. 



3. Catalogue cross association 

The starting point for the cross-association process is this set 
of individual SPIRE band catalogues. For the current work the 
emphasis is on the creation of a robust, un-confused sample of 
sources that has the highest probability for its colours to origi- 
nate from single unblended galaxies. We only consider the cen- 
tral regions of the maps, that have full homogeneous cover- 
age by all scans. To protect against spurious sources, we com- 
pare fluxes separately derived from two independent half-maps. 
The ratio of the two flux estimates separates well into 3 dis- 
tributions. Spurious sources are removed by excluding ratios 
above 5 and below 1/5. For this work we have constructed a 
500 jt/m band flux limited selection. It is justified in three ways: 
i) The stronger negative /T-correction means selection in this 
band favours higher redshift galaxies; ii) This is a relatively 
new band as yet only ex plored by much shallower BLAST sur- 
veys (iDevhn et al.ll2009h : iii) About ~ 90% of 500 jum selected 
sources are also detected in the other SPIRE bands. We require a 
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Fig. 1. Examples of single 500 yum source detections {top and 
bottom panel, right) with multiple counterparts at 350 ;um and 
250 //m {top panel, middle and left), and single 350 and 250 jum 
counterparts {bottom panel, middle and left). Note that the up- 
per 500 jum source appears already double to the eye. The red 
marked source position in each image is the one determined by 
the algorithm at 500 jum. Each image measures 3' on each side. 
The yellow circles indicate the FWHM of the beams. 



signal-to-total-noise (S/N) ratio of more than 3 in the 500 //m fil- 
ter. The formal average flux uncertainties at 500 yum derived from 
the source extractor results are only 1.1, 2.1, 4.4, and 11.6 mJy. 
We consider these to be instrumental noise, based on their ratio, 
consistent with the 30, 7, 2, and 1 repetitions executed on the 
four fields, respectively, and the 1 cr confusion noise at 500 /vm of 
6.8+0.4 reported bv lNguven et alJ (1201 0). Thus the uncertainties 
are completely dominated by extragalactic confusion for the first 
two fields and increased for FLS and Lockman-SWIRE. We cal- 
culate the following effective 3cr average flux limits in our fields 
from 3 times the average total eiTor of all sources with S/N< 4: 
23.4, 26.4, 32.0, and 46.4 mJy for GOODS-N, Lockman-North, 
FLS, and Lockman-SWIRE respectively. The selection leaves 
48, 61, 608, and 824 sources at 500 //m in the 4 bands respec- 
tively. This conservative threshold also minimises the impact of 
flux boosting on the derived colours of the sources. 

To further de-blend and cross-match, first, all 500 jum 
sources without another 500 yum source within an 18" radius are 
selected. This radius was chosen to be similar to the beam size 
at 500 jjva. Then for these remaining sources, the same 18" ra- 
dius is checked in the other two bands. Sources with more than 
one source in a different band are discarded immediately. In case 
only one source is found in the other band, it must be within a 
radius of 8" in order to be accepted as a cross identification, 
otherwise the source is considered blended and discarded. This 
radius was chosen to include 3-cr of the telescope pointing error 
and the estimated PRE fit error of 6" each. We end up with a list 
of potentially uncontaminated 500 /im sources that is then cross 
matched with the lists of the other two bands with a match radius 
of 8". In Figure [T]the dangers of simple naive associations are 
emphasised, where a cluster of sources shows up as being de- 
tected as single at 500 //m by the point source extractor, but re- 
vealing multiple counterparts at 250 yum. This sample is largely 
free from contamination and should have reliable fluxes originat- 
ing from just one source, accurate at a 30 % level or better. The 
final matched source numbers for the four fields respectively are 
21,38,242, and 244. 



4. Analysis 

In Figure |2] we plot the 3-dimensional SPIRE flux-flux-flux pa- 
rameter space for our band merged catalogue. The fluxes are 
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Fig. 2. The 3-dimensional flux parameter space for our un- 
blended band-merged catalogues in the SPIRE 250, 350 & 
500 yum bands. The sources in GOODS-N appear in red, 
Lockman-North in blue, FLS in grey, and Lockman-SWIRE in 
green. Both diagrams show the same 3D plot from two different 
aspect angles, the left one from within a plane fitted through the 
data, the other from a perpendicular direction. 



grouped around a relatively flat and thin surface in the 250 yum, 
350 ;um, 500 ;um parameter space. The same even thinner surface 
is seen in similar plots of mock catalogue data that is discussed 
later. Thus, although we have flux data in three SPIRE bands, 
in principle only two parameters are needed to describe the in- 
formation. This degeneracy follows from the fact that the spec- 
tral energy distributions (SEDs) in the submm, which SPIRE 
observes, are dominated by dust emission that have very simi- 
lar shapes and result in fairly well defined flux ratios. Thus, the 
main parameters determining the three SPIRE fluxes, are rather 
wavelength of the emission peak and luminosity. 



4.1. Colour - colour parameter space and model comparison 

In the Figure |3] the Siso/Ssso - S350/S500 colour-colour diagrams 
for SPIRE sources in the SDP survey fields are shown with the 
co lour tracks from the contemporary galaxy e volution r nodels 
of iPearson et al.' ("200^, iDale & Helour(l2002h . IXu et al.l (l200lh 



flPe 



&|Lagache et al. (2003!) over-plotted on individual panels. The 
redshift of the tracks is shown in colour and ranges from to 4. 
In general, all the models are consistent with the obtained 
SPIRE colours, except from individual su btleties of the models . 
However, especially the SED templates of iPearson et al.l (l2007h 
and Lagache et al. (2003) lack diversity in dust illumination con- 
ditions to cover the spread of colours sufficiently. All models 
agree: the colours imply that the SPIRE population is not lo- 
cal but rather the bulk hes at reds hifts between 1 and 3 .5. We 
note that the IXu et al.l (l2001b and iDale & Heloul (l2002h mod- 
els tend to place the population at som ewhat lower redshif t 
tha n the other two. Thi s implies that the IPearson et al.l (l2007h 
and lLagache et al.l ( 12003 ) model SEDs contain generally warmer 
dust, which is confirmed by plots of the emission maxima of the 
SEDs. 



four crosses on the left are in the same vertical order as the sym- 
bols in the legend and represent the averaged uncertainties in 
the four fields. Different tick marks show instrumental and total 
components. The four vertical lines indicate from left to right, 
the eff'ective flux limits of GOODS-N, Lockman-North, FLS, 
and Lockman-SWIRE respectively. In both panels of Figure |4] 
the observed data are compared to mock catalogues of 1 deg^ 
on the sky by Pearson et al. (2007) (left) and Xu et al. (200ll) 
(right) that were cut below the eff'ective flux limit of GOODS- 
N. Again the most notable difference is the larger spread of the 
IXu et al.l (2001) colours due to a larger number and diversity of 
SED models. In both models the bulk of objects are Starburst 
galaxies, LIRGs and ULIRGs, that are grouped around a colour 
of S250/S350 ~1-1- 

The high-redshift sources populate a specific area of the 
colour plane in both models, although the redshift distributions 
are different. In the IPearson et al.l (|2007|) model the highest red- 
shift objects, z>3 occupy the parameter space corresp onding to 
S2.50/S 350 colours <L0 with 8500 <40mJy, while the IXu et al] 
(I2OOI I) model locates the z>3 region rather at S250/S350 <0.8 and 
the same flux cutoff, but with fewer objects and mixed with many 
low redshift SEDs. Similar cuts can be made for z>2 sources. 
Lower redshift sources may also be excluded by virtue of their 
higher S250/S350 colours. 

The SPIRE data generally overlap fine at S5oo<60, except for 
colours of S250/S350 < 0.8. Especially the Pearson et al. (2003) 
model shows no objects below this Umit , while the same re- 
gion is sparsely populated by the IXu et al.l (1200 lb model SEDs. 
Looking at the model types, it turns out that those are mainly 
AGN, which are missing entirely in the ' Pearson et al.l (l2007l) 
SED catalogue. Neither model covers sufficiently the increas- 
ingly redder colours in this region that SPIRE observes towards 
500 jum fluxes above 50 mJy. A comparison with another mock 
catalogue by Valiante et al. (2009) shows the same lack of red 
sources. 

A considerable fraction of submm bright s ources are ex- 
pected to be lensed by foreground galaxies dNegrello et al.l 
L2007). Since lensing magnification is wavelength independent, 
such lensed sources appear in their intrinsic positions in the 
colour-colour diagram, but their locations in the colour-flux 
plane would be off'set to brighter fluxes (towards the right side of 
Fig. ID along t he x-axis), whil e keeping col ours the same. Using 
the m odels of Negrell o et al.l (1200 7) and Negrello priv. comrnj 
(I2OI O) we estimate that of all objects with fluxes S500 > 100 mJy 
and redshift 2<z<3, almost all are lensed. For our 19.6 deg^ total 
sky area, that would be ^2-3 out of the 24 bright 500 yum sources 
we have identified, although this is a lower limit, as our selection 
procedure may bias slightly against clustered and as such poten- 
tially lensed sources. 

For now we conclude that we see a population of red bright 
objects that may consist mostly of colder SEDs but with a frac- 
tion of di stant lensed ones. Inclusion o f other wavelengths as 
shown by iRowan-Robinson et al.l (1201 Ol) will be needed for fur- 
ther interpretation. 



4.2. Colour - flux parameter space 

To overcome the apparent degeneracy in the SPIRE colours 
(Figure |2]), we plot colour-flux distributions. In Figure |4] we 
show the S250/S350 colour versus 500 /im flux density distribu- 
tions for the SPIRE sources. A few sources at S5()o>100 mJy are 
not shown to improve visibility. The symbols indicate the differ- 
ent fields according to the legend in the upper left comer. The 
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Fig. 3. Ssso/Ssoo - S250/S350 colour-colour plots f or the SP I RE so urces. Over plotted are the colour tracks from the g alaxy evolution 
models of top-Ze/f Pearson et al. (2007), top-risht \Xu et al.l(l2001h . bottom-le ffiLa^ache et al.l(l2003l) . and bottom-risht \Dale & Heloul 
(l2002ft . The redshift in the tracks is colour coded and runs from to 4. The black symbols represent all unblended SPIRE sources 
according to the legend. 
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Fig. 4. Meas ured S^.sn/S^'in colour 500 fim flux distributions fo r the SPIRE sources (black symbols) in comparison with mock 
catalogues of iPearson et al.l (l2007h to the left and lXu et al.l (1200 Ih on the right. The large error crosses on the left represent average 
1-cr total uncertainties dominated by extragalactic confusion and the smaller tick marks show instrumental noise only, which is 
negligible for GOODS-N and Lockman-N. The four vertical lines indicate from left to right, the effective flux limits of GOODS-N, 
Lockman-North, FLS, and Lockman-SWIRE respectively. 
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